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The global behavior of scalar field cosmological models 
with very hard potential walls is investigated via the simple 
example of an exponentially steep potential well. It is found 
that the solutions exhibit a non-trivial oscillatory behavior 
in the approach to an initial space-time singularity. This be- 
havior can be interpreted as being due to the inability of the 
scalar field to negotiate the walls of the steep potential well. 



I. INTRODUCTION 

According to the inflationary universe scenario 
the matter content in the very early universe can be mod- 
eled by a scalar field with a non-negative self interaction 
V{(f>). Inflation may be understood physically in terms 
of the stress energy tensor of the scalar field which is 
equivalent to a perfect fluid with energy density p and 
pressure p given by 
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When V becomes large compared to the gradient of 
the scalar field, the pressure becomes large and nega- 
tive and the fluid satisfies the approximate equation of 
state p — — p leading to a rapid exponential expansion 
of space-time. This rapid expansion is associated with 
the fact that gravity becomes repulsive in the presence 
of large negative pressures. Such equations of state have 
traditionally been forbidden in general relativistic matter 
fields by the strong energy condition which, for a perfect 
fluid, is equivalent to the following familiar constraint on 
the pressure and energy Q: 

p + 3p>0. 

Thus, for a scalar field to obey the strong energy condi- 
tion we would require that 

-d^<f>d^^ - v{(j>) > 0. 

The fact that scalar fields naturally violate the above 
inequality means that the singularity theorems Q can 
not be strictly applied when such fields are present in 
the early universe. Indeed, it is relatively simple to con- 
struct exact singularity free scalar field cosmologies (eg 
1^). Other exact scalar field solutions are known which 
possess singularities, but have no particle horizons 
indicating that scalar field cosmologies can exhibit a di- 
verse range of behavior. 



It turns out however, that all of the known singu- 
larity free and horizon free solutions are unstable and 
it can be demonstrated that, provided the potential di- 
verges slower than exp(\/60) as (f> ^ ±oo, almost all spa- 
tially flat Friedmann Robertson Walker (FRW) cosmolo- 
gies have initial singularities. Furthermore, they have 
asymptotic equation of state p — p and may be asymp- 
totically approximated by the general solution for the 
massless scalar field . 

However, since we have no direct emperical data with 
which to determine the detailed nature of physical fields 
at energy densities comparable with those of the early 
universe, we should not rule out the possibility that scalar 
fields present in the early universe might have much 
steeper potentials. 

It is therefore natural to ask whether singularity free 
or horizon free cosmologies can arise as a result of very 
steep potentials. In particular, we might expect that a 
very steep potential well could inhibit the divergence of 
the scalar field, thereby slowing down the gravitational 
collapse and resulting in singularity or particle horizon 
avoidance. 

Indeed, it can easily be demonstrated that for very 
steep potentials V can not be neglected asymptotically 
in the past as can be done with less steep potentials. This 
can be seen as follows: 

For simplicity, we confine attention to the FRW line 
element 
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We define the expansion K{t) as the function 



K 
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where the "dot" indicates differentiation with respect to 
time t. K may be interpreted physically as the rate of 
expansion of the spatial volume element v = a'^. 

The evolution equations governing the interaction of a 
general scalar field with the above metric are M : 



4, = -K^-V'{cj)) 
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The general solution for the massless scalar field is ob- 
tained by setting 1/ = in the above and solving to 
obtain. 
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Let us now use as an approximation to the solution 
for arbitrary V in the hmit t ^ 0. Substituting into 
(H) we obtain 
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where, /i indicates terms of higher order in t. For this 
expression to be consistent the foUowing must be true: 



UmtV ±A/|hi 
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I.e. 
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In other words, in order for the massless scalar field 
to represent a consistent asymptotic approximation of 
a scalar field model it is necessary that the potential di- 
verge slower than e^^'"^'. Exponential potentials there- 
fore seem to accommodate a kind of transition between 
soft wall potentials and some other qualitative regime. 

In this paper we examine the behavior of very steep 
potentials by means of a simple and interesting example. 



II. THE MODEL AND DYNAMICAL 
EQUATIONS. 

The simplest example of a potential which is exponen- 
tially steep at both infinity and negative infinity is an 
exponential well of the form. 

V{(p) = ae^-^ + be~'"t'. 

where a, 5, A and fi are positive constants. In the subse- 
quent analysis we will confine attention to the case where 
a = b = 1 and A = /i so that 



(7) 



By increasing A we shall be able to investigate how the 
qualitative behavior of the system changes as the poten- 
tial well becomes increasingly steep. The analysis for the 
more general potential is similar and the conclusions are 
essentially the same. 

We are interested only in expanding solutions of (|^) so 
we may replace K and 4> with the new set of coordinates 



X 
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and introduce a new time coordinate 
T = \nv{t) + To 



(8) 
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where w = a'^ is the spatial volume element and tq is some 
constant, r is well defined since v is strictly increasing. 
Furthermore, u goes to zero if and only if K goes to 
infinity so (ph implies that i^T— > ooasr^ — oo. 
Differentiating ffl) we find (recalling that K = v/v) 



d d 

-r = K—. 

dt dr 
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In terms of these coordinates the field equations (Q) 
with the potential (^ may be written as a dynamical 
system: 



dx 
d^ 
dy 

dr 
d(j> 
d7 
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where a = ^ . The "constraint" equation (|5|) becomes 

2/2-H3a;2(e^"'^-he-^"'^) = 1. (12) 
We can simplify the equations by defining the variables 



p = V3e-^°"''x 



q = y- 



then the constraint equation may be written 



(13) 
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Substituting (|13D and (|lj) into (|l|) we obtain 

dp 2 
-_ = -apq + pq 



dq 3 
d7 = '^ 



aq'^ — q — a + 2ap^ . 



(15) 



These equations constitute a 2-dimensional dynamical 
system on the p-q plane. We may define the physi- 
cal phase space 17, according to (p^, as the set fl — 
{{p,q) : + q^ < l,p > 0}. In other words, all physi- 
cal trajectories lie on the interior of the unit disc to the 
right of the q axis. The unphysical boundary dfl is a 
closed curve consisting of the union of the smooth arc 
dfli — {{p,q) : p^ + q^ = l,p > 0} and the line seg- 
ment — {{p,q) ■ P — ^ 1}- The asymmetric 
appearance of dfl is a consequence of the positive ex- 
ponential term in the definition of the coordinate p and 
is not a physical property of the system itself. In fact it 
should be pointed out that dQi maps onto 8^2 under the 
transformation {(f>,(f>) ~{<j),(j>). dfl corresponds to the 
infinity of the expansion K. In order to see this define 
the function 



i/(p,(?)=p2(i„p2_^2) 



(16) 



Clearly H is strictly positive everywhere on O and van- 
ishes identically on dfl. From (|l3) and (Oh we see that 
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Thus do, is just the set of ah points for which x — 
which is by definition the infinity of K . Evaluating the 
directional derivative of H along the flow using ( [l5| ) we 
find 
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which is non-negative everywhere on the interior of . 
In fact the derivative of H with respect to r is strictly 
positive everywhere on the interior of except where 

dH 

q = Q. Observe also that — — = on dVt indicating that 

dr 

the boundary is an invariant manifold (tangent to the 
fiow). No physical trajectories can therefore cross dVl 
into the unphysical domain beyond. Since f2 is compact 
all trajectories must possess an past-limit set which is in- 
variant under the flow. Since H is monotonic, any limit 
dH 

point must have — — = 0. Therefore, all limit sets must 

be subsets of either dVl or the line q ~ Q. From ( |15| ) the 
only invariant subset of g = is the equilibrium point 
■pd = (-^j 0)0. However this point is a local maximum of 
H as is easily verified by evaluating its gradient. Since H 
is monotonic increasing, no solutions can be past asymp- 
totic to pd other than the steady state solution on 
itself. It follows that the past-limit sets of all other solu- 
tions lie on the boundary dVl. 

By the time reverse of the above argument it is clear 
that all solutions, with the exception of those unphysical 
solutions lying on 917 are future asymptotic to p^. 

The future asymptotic set pd represents a vacuum de 
Sitter space-time with constant expansion K — \/6 and 
(j) identically zero. This is consistent with what we would 
expect for a scalar field cosmology with non-zero vacuum 
energy. 

To summarize, it has been established that (almost) 
all solutions originate near dU, and subsequently evolve 
towards the global attractor p^ (which represents de Sit- 
ter space) as i cxD. Let us now examine the behavior 
of the system on and near in more detail. 



III. THE BEHAVIOR CLOSE TO dVl. 



There are a maximum 4 equilibrium points of (15) ly- 
ing on aO. These are (a/(1 - ^2), a), (0,1), (0,-1) and 
(0, —a), which we label pi, p+, p_ and P2 respectively. 
Observe that pi and p2 only exist as distinct equilibrium 
points on 917 when a < 1. For values of a greater than 
or equal to one, p± are the only equilibrium points. In 



order to investigate the behavior close to p± equations 
(15) and (16) can be linearized about these points. The 
linearized system is: 



dp 
dq 



(1 T a)p 
= 2(l±a)(9Tl) 
The solution to the linear system is 



(19) 
(20) 



p^poe^^^"^^ g = ±1 ^ ^oe2(i±")" (21) 

where pq and Sq are positive constants. 



A. The Flow For a < 1. 

When a < 1 both exponential terms have positive co- 
efficients indicating that p± are sources of the linear sys- 
tem and therefore, by the Hartman-Grobman Theorem 
§,|, of the non-linear system also. Since the unphysical 
solutions on dfl move away from p± they must asymptot- 
ically approach pi and P2 (in the forwards time sense). 
However all physical solutions (ie, those on the interior 
of 17) asymptotically approach pd in the future, pi and 
P2 must therefore be unstable saddles. This can be seen 
more clearly from inspection of the geometry of the flow 
as illustrated in Fig. |l|. 

( An alternative way to verify that pi and P2 are saddles 
is by linearizing ( [l5|) about these points but this exercise 
shall be left to the reader who remains unconvinced). 




FIG. 1. Sketch of f2 showing all possible equilibrium points, 
direction of the flow on the boundary, and some typical tra- 
jectories. Trajectories on the boundary are future asymptotic 
to pi and p2 but all trajectories on the interior approach p^. 



^In order to avoid confusion points on the (p, g) plane will 
be labeled in bold print. 
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With the exception of the 2 solutions past asymptotic 
to pi and P2 respectively, and the steady state solution 
lying on p^, all physical trajectories must asymptotically 
approach either p+ or p_ as r — > — oo, as indicated in 
Fig. |ll 

Using the fact that 1 — decays exponentially to zero 
as r — > —oo we may integrate (ns) to obtain a first order 
expression for H: 



H = Hoe 



4r 



and hence 



X = XqC . 

Using ( |lo|) we obtain a first order expression for t: 

t = XqC'^ . 

Observe that t^Oasr— >— oo indicating that p± cor- 
respond to space-time singularities. Using the definitions 
of X and q we thus obtain the asymptotic solution for K 
and (j) in the neighbourhood of < = 0: 



K 



± 
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(22) 



terms, the gravitational expansion is unable to dominate 
the scalar field self interaction in the initial expansion 
phase of the universe. For expediency, we shall assume 
in what follows that a > 1. The special case a = 1 differs 
in the details of analysis, but the qualitative features of 
the solutions turn out to be essentially the same. 

The only equilibrium points possessed by the system 
when a > 1 are p+, p_ and the sink Poo . Inspection of 
( pl| ) reveals that p± become hyperbolic saddles in this 
regime. 

The only solution past asymptotic to p_ is the unphys- 
ical trajectory 71 which emerges from p_ at r = —00 and 
proceeds anticlockwise along d^li (the unit circle), reach- 
ing p+ at T = 00. Similarly, the only solution which 
originates at p+ is the unphysical trajectory 72 which 
emerges from p_(- at r = — 00 and proceeds along 8^2 
(the g-axis), approaching p2 asymptotically as r — s- cx). 
The union of 71 , 72 is a closed heteroclinic cycle 7l on dfl 
which forms a limit cycle for trajectories on the interior of 
fl. That is, as ?7 ^ —00 a typical trajectory will approach 
dQ, spiraling clockwise (in the reverse time sense) an in- 
finite number of times. The solution will asymptotically 
approach the non-physical solutions 71 and 72 but unlike 
these will always avoid the equilibrium points p+ and p_ 
and will continue to spiral around dQ ad infinitum. Fig 



and upon integration of 1 
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This is, of course, the general solution for the massless 
scalar field. We thus conclude that when a < 1, ie A < 
VQ, the potential V is not dynamically significant near 
the singularity. This is precisely the behavior we would 
have expected from the simple calculation in Section 1 for 
potentials that go to infinity slower than exp{^/6(l>). In 
fact the topological structure of solution space is typical 
of inflationary models. All s solutions emerge from the 
sources p± except for the seperatrices Ii or I2 which 
we define as the unique solutions past asymptotic to the 
saddles pi and p2 respectively. We may infer that these 
solutions characterize the inflation in the system. (It 
can be easily verified that the behavior close to pi and 

t2 is approximated by the exact power law cosmologies 
|. However this behavior only corrosponds to infiation 
when < i since this condition ensures that the square 
of the q coordinate at p+ and p_ is less than i which is 
necessary and sufficient for violation of the strong energy 
condition.) 




FIG. 2. Sketch of Q for a > 1 showing the solutions on the 
boundary and a typical solution on the interior. The points 
p± are saddles and therefore unstable, but solutions must 
approach the boundary so the union of solutions on dil forms 
a limit cycle. 



B. The Flow For a > 1 

Let us now consider what happens when a > 1. As 
we demonstrated in Section 1, the potential in this case 
becomes too steep for the solution with F = to be 
consistent as an asymptotic solution. In more physical 



In order to interpret these results recall firstly that, 
from ( p^ ) and (|l^), H monotonically decreases to as 
T —^ —00 and hence the expansion K diverges monotoni- 
cally to infinity. The asymptotic periodic behavior of the 
trajectories must therefore represent oscillations of the 
scalar field (j). Inspection of (|lj) indicates that cj) = —00 
on the semi-circle dfli. It follows from symmetry that 
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= cx) on 80,2- As the gravitational expansion diverges 
to infinity ( and correspondingly the scale parameter a 
tends to 0) the scalar field oscillates about (j> — with 
the amplitude of oscillation increasing for each successive 
cycle, asymptotically diverging to infinity. 

No asymptotic equation of state is satisfied by the 
scalar field. In the "corners" of fl, close to the points 
p± , the equation of state is approximately p ~ p since 
02 ~ which imphes, via (||) that 0^ >> V{<j)). This 
corresponds to the period where the scalar field is close 
to the bottom of the potential well. 

Close to the g-axis, on the other hand, the equation 
of state is approximately p — —p since </) ~ whilst 
Vif/)) ^ 2. This corresponds to the scalar field reaching 
the top of its roll before accelerating down the potential 
well again. The asymptotic behavior of the space-time 
may therefore be characterized by an infinite sequence of 
"stiff" phases, punctuated by de Sitter phases. During 
the de Sitter phase the matter violates the strong energy 
condition. 



IV. THE EXISTENCE OF A SINGULARITY AND 
PARTICLE HORIZON 

Since an infinite number of physical cycles of the scalar 
field occur for the case a > 1, it might reasonably be 
expected that an infinite interval of proper time must 
also elapse since the oscillations of the scalar field could 
be used as a natural physical clock which would measure 
an infinite time interval to the past of any space-time 
point. If this were the case then the space-time would 
be non-singular for all physical solutions. It turns out 
however, that the oscillations pile up on each other and 
a space-time singularity is indeed reached after a finite 
time interval. 

In order to show this we must estimate the proper time 
that elapses to the past of some arbitrary initial point 
Po, on the trajectory ippg^r). Since all trajectories are 
asymptotically tangent to the boundary dVl it will be 
sufficient to restrict attention to initial points lying in 
the set TiQ — {{p, q) ■ q — I — €, < p < e} where e > 
may be chosen arbitrarily small. 

Before proving that space-time singularities and parti- 
cle horizons exist for all solutions it will be convenient to 
prove the following lemma; 

Lemma 1 For all < n < I there exist positive numbers 
e andpm such that if Xq = (po, 1 — e) is in Sq andpo < pm 
then 



H{^p„iT)) < Hoe 



(24) 



to the past of po , where ip^o (t) is the unique trajectory 
of (|7^ passing through po with ipp^ (0) = po . 



Proof: Integrating (|1^) backwards in time from r = to 
some earlier time Tf we obtain 



iJ(r/) = iJoexp -4/ q^dr 
Let El = {{p, q) e ^ : p = €, 1 - g < e}. 



(25) 




FIG. 3. The box Q,t . On each successive cycle in the re- 
verse time direction the trajectory intersects the box for a 
finite period of time coming progressively closer to the bound- 
ary of f2 (the arrows indicate the direction of the flow in the 
forward time sense ) . 

As can be seen from Fig. ^, the set Ei -|- Eq encloses 
a box, f2+, of area ~ in around p+. Another box, 
, can similarly be constructed around p_ by defining 
the sets E2 = {(p, q) € ^ : p = t 1 + q < e} and E3 = 
{(p,g) :q = -l + e, 0<p<e}. 

Let / be the time interval [t/,0], then for a given tra- 
jectory V'po we may write / = /c U /{, where /c = {r S / : 
^Po(^) G ^f} and /fc = {r e / : V^p„(t) ^ fif }. We thus 
have 



q^dT = / q'^dr + / q'^dr 
> q^dr. 



(26) 
(27) 



From the definition of Ir 



q^dr = (1 - 0{e))dT. 



Fix n and let m be any number satisfying n < m < 1. 
Then for e sufficiently small we have; 



q dr > m / dr. 



(28) 



In order to complete the proof we must show that by 
choosing po sufficiently small, the ratio of Jj dr to Jj dr 
may be made arbitrarily close to 1. 
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According to (jlj) p evolves according to the equation 

d 



dr 



Inp - 



-aq + q 



<l + a 



since \q\ < 1. Therefore if T2 < ti we have 
In^ < (l + a)(Ti -T2), 

P2 



(29) 



(30) 



where pi — p{ti),P2 = p{t2)- The parameter time At 
needed for a point po = (po, 1 — e) G So to flow through 
and reach Ei must therefore satisfy 



lArl > 



1 



1 



hr - 



(31) 



The modulus sign is necessary because we are tracing 
the trajectory backwards in time so Ar is negative. As 
Po — > , Ar — > — 00. Thus Ar may be made arbitrarily 
large by choosing po sufficiently small. 

Since "i/^pp is past asymptotic to the limit cycle 7;, the 
intersection of the past orbit with the line q — 1 — e 
(which contains So) must contain an infinite number of 
points in addition to po itself. If pi is any such point 
then substituting g = 1 — e into ([l^) and using the mono- 
tonicity of H dlq) we must have pi < po, where pi and 
Po are the p-coordinate values of pi and po respectively. 
In other words ^pp^ must intersect E^o again after 1 com- 
plete cycle of dfl and the point of intersection pi must 
have a p-coordinate value pi which is smaller than po. It 
follows that V'po passes through the box ri+ on each suc- 
cessive cycle and the time interval Ar+ to traverse 
on the rt"* cycle obeys the inequality; 



a + 1 \Pn 

> In — 

a + 1 \po 



(32) 
(33) 



Similarly, defining Ar~ to be the time taken for ipp^^to 
traverse the box on the n*'' cycle we find from (^) 
using an identical argument to that above that 



IAt: 



> 



— In — 

a+1 \Pn 



(34) 



where p„ is the p coordinate of the intersection of ippg 
with S3 on the n*'' cycle (by the n*'' cycle I mean, pre- 
cisely, one complete circuit from (p„ — 1,1 — e) e Sq 
to (p„, 1 — e) € So). Recalling the definition of H 
and using the fact that q^ takes the same value on 
So and S3 we have for any point on So or S3 that 
= (2e— e^)p^-l-O(ep^). Since i7 is monotonic it follows 
that for e chosen sufficiently small p„ < po- Therefore, 



|Ar-| > 



a+l \po 



(35) 



Let us now consider the parameter time Ar taken for a 
point on Si to reach S2 under (|l5|), close to the semi- 
circular boundary dili. That is , the time taken to flow 
backwards in time from the top box to the bottom 
box (Fig. I). 




FIG. 4. The parts of the solution flowing between the boxes 
and $17 approach the restrictions of 71 and 72 to finite 
time intervals. 

By continuity, Ar approaches the (negative) param- 
eter time interval which the asymptotic solution 71 
takes to map the point (e, Vl — e^) G Si to the point 
(e, — Vl — G S2. This interval is finite. 

Therefore, the time interval to go from Si to S2 must 
approach a finite (negative) limit. Its modulus must 
therefore possess a finite upper bound ri . 

Let Ar^ be the time taken for tp^o to go from Si to 
S2 on the n*'' cycle. Then we have 



|Ari|<ri 



(36) 



Similarly, if Ar^ is the time taken for tjjj;g to go backwards 
in time from S3 to So on the n'^ cycle we have 



|Ar^|<r2 



(37) 



For some finite number r2. 

Now, Since is incident on So and hence, initially 
fiows through flf , it follows from the inequalities js^), 
(H), (H) and (^) that 

Jj^dr re (a -1-1) \po 
where r^ = max(ri,T2). Choose 



Po = e exp 



(f -1) 
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then we have 





dr 


{•' 








771 




dr. 




L 









(38) 



Combining (|38|), (|2§) and (g5|) gives 
H{Tf) < Hoe^^'^f 
which is the required result. □ 

We say a scalar field cosmology is non-trivial if there 
exists some space-time point for which is non-zero. 

Theorem 2 // (g^^j/ , 0) is a non-trivial scalar field cos- 
mology with potential ^ and if g^^ is spatially flat and 
isotropic, then 5^^^ possesses an initial space-time singu- 
larity. 

Proof: 

By Lemm a u\ there exists e and pm such that all tra- 
jectories of dlq) incident on Eq with pq < pm satisfy 

H{t) < Hoe^^ 

for all T < 0. It will be sufficient to show that these tra- 
jectories reach the boundary in finite proper time. The 
proper time taken to reach the boundary is given by 



At: 



xdr 



^jl^dT 



< Xo 

= 2x0 

Since this is finite, we have the result. □ 



(39) 

(40) 
(41) 



Theorem 3 // {g^iuT4') is a non-trivial scalar field cos- 
mology with potential ^ and if g^y is spatially flat and 
isotropic, then particle horizons exist for all isotropic ob- 
servers (observers whose worldlines are tangent to the 
timelike Killing vector field U^). 

Proof: A particle horizon exists for an isotropic observer 
at time t if the integral 



I = 



1 



dt 



exists and is finite. By the definition of r; 

a = 

and 



dt 
d^ 



Therefore, 



e xdr. 



By Lemma y there exists e and Pm such that all trajec- 
tories of (nS ) incident on Sq with po < Pm satisfy 



H{t) < Hoei^ 



(42) 



for all r < 0. If Hs finite at r = it will be finite for all 
T. It will therefore be sufficient to show that trajectories 
for which (^2|) holds possess a horizon at r = 0. Using 
(E^) we have 



/ < Xq 



e^dr 



(43) 

= 6x0 (44) 
Since this is always finite we have the result. □ 



Note also that as xq 0, I must also approach indi- 
cating that the horizon length shrinks to zero as t — > 0. 

The physical meaning of the Lemma |l| becomes clearer 
when we recall that r = Inz; and, by and the defini- 
tion oi X H — QK^'^. Lemma ^ may thus be interpreted 
as saying that for all < 71 < 1 there exists some A> Q 
such that 

K > Aw"". 

Recalling that K = ii/v it thus follows that for all 777 > 1 
there exists vq such that 



V{t) > VQt" 



(45) 



Thus, in the neighbourhood of the singularity the volume 
element expands faster than any power law with power 
greater than one. In order to avoid a particle horizon it 
must expand slower than t^. 

What about the case 777 = 1? Consider the function 
g-i'^jj Taking the derivative of this function with re- 
spect to T, using 

d 



(18) 



dr 



e-^^H = -4e-*^i/(l - g^) 



Thus for r < we have 



.0 

In (e-'^^H) =c + 4: (1 - q^)dT 

where c is a constant. The integral on the right hand side 
tends to infinity as r — > — 00 since each solution spends a 
finitely large amount of r-time with, say, < e on each 
cycle (of which there are an infinite number). 



7 



Thus, 



lim 



''H = oo. 



This translates to a corresponding hmit for v and t, as 
above: 



hm t^^v = 0. 



(46) 



Comparing this expression with ( |4^ ) we see that v ex- 
pands slower than t but faster than any power law 
with p > 1. In this respect the behavior of the gravita- 
tional field near the singularity is quite subtle and un- 
usual since it can not be adequately modeled by a power 
law. Note also that it is clearly not admissible to neglect 
the dynamical effect of the potential when considering 
the gravitational field near the singularity. 



V. CONCLUSIONS 

The above example indicates that non-trivial asymp- 
totic behavior can emerge, even in very simple models, 
as a result of very steep self interaction potentials. It is 
of interest that both space-time singularity and particle 
horizons exist even though the strong energy condition 
is violated by typical solutions during all periods of their 
evolution, including asymptotically close to the singular- 
ity. 

It is suggested by the author that the qualitative fea- 
tures exhibited above, including the oscillatory behavior 
and existence of a singularity and and particle horizons, 
are characteristic of FRW scalar field models which have 
very steep potential wells. Preliminary investigations of 
the steeper than exponential potential = e^'^ and 

the "hard wall" potential V{4>) = 1/(A - cj}"^) on the do- 
main 0^ < A reveal no significant departure from the 
qualitative behavior of the exponential potential well. 

The asymptotic oscillatory behavior of the scalar field 
is reminiscent of the oscillatory behavior of the com- 
ponents of the shear tensor displayed by the Belinskii, 
Khalatnikov, Lifschitz perfect fluid cosmologies 
This behavior is associated with the existence of dynam- 
ical chaos in the solution space of these cosmologies . 
The apparent similarity suggests that it might be par- 
ticularly interesting to investigate the dynamics of ex- 
ponential potential well models in slightly more compli- 
cated space-times such as Bianchi type IX. It seems not 
inconcievable that the coupling of scalar field and shear 
degrees of freedom could increase the degree of mixing, 
thereby resulting in partical horizon avoidance, particu- 
larly since, unlike the BKL solutions, the expansion of 
scalar field cosmologies is not constrained by the strong 
energy condition. 

Another interesting feature of the solutions is that 
there exists two natural time scales, namely the affine pa- 
rameter time t and the period of oscillation of the scalar 



field, which we might choose to characterise by some co- 
ordinate 77. The timelike geodesies are incomplete with 
respect to t but complete with respect rj. 

By convention we say that a space-time singularity ex- 
ists if geodesies are incomplete with respect to their affine 
parameter since this describes a situation where a freely 
falling observer would reach the edge of space-time in fi- 
nite proper time. However, the physical interpretation 
of t as the proper time is not really meaningful on a 
neighbourhood of a boundary point of space-time since 
no normal coordinates can be constructed at such a point 
(normal coordinates can be constructed at a point arbi- 
trarily close to a singularity but may never be extended 
to the singularity itself). 

On the other hand, rj can be interpreted as the most 
natural time scale associated with the matter content 
of the universe (including all clocks and astronauts since 
scalar matter is the only matter entering into this model) 
and therefore might be a more meaningful measure of 
the time taken to reach the singularity. One could imag- 
ine that if a more realistic cosmological model could be 
shown to display a similar oscillatory behavior for its 
physical fields then one would be led to an interpretation 
whereby a singularity exists but can never be reached in 
a finite amount of time by any physical object. 
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